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1. Introduction 

In this paper we characterize the ground state energy of hydrogenic atoms in 
magnetic fields. We deal with fields of large strength in the Dirac operator frame- 
work, far away from the perturbative regime. 

To compute eigenvalues of Dirac operators, the usual min-max principle does 
not apply. More sophisticated versions of this principle have been established 
over the last few years, see lllJlUl!3| . These techniques are powerful enough to 
provide accurate and efficient algorithms for calculating eigenvalues of Dirac 
operators 7,6 . In this paper we demonstrate that they are also flexible enough 
to cover the case with a magnetic field and provide reasonable results for a highly 
non-perturbative problem, when paired with the right physical insight. 

The Dirac operator for a hydrogenic atom in the presence of a constant mag- 
netic field B in the xa-direction is given by 



Hb — 1 — r with Hb ■= a 

\x\ 



11. 

- V-t- -B(-X2,xi,0) 
I 2 



(1) 



where u = Za < 1, Z is the nuclear charge number. The Sommerfcld fine- 
structure constant is a ~ 1/137.037. The energy is measured in units of mc^, 
i.e., the rest energy of the electron, the length in units oih/mc, i.e., the Compton 
wavelength divided by 2 tt, and the magnetic field strength B is measured in units 

of ™^|^ « 4.4 X 10 Tesla. Here m is the mass of the electron, c the speed of light, 

q the charge of the electron (measured in Coulomb) and h is Planck's constant 
divided by 27r. It is worth recalling the the earth's magnetic field is of the order 
of 1 Gauss and 1 Tesla is 10"* Gauss. 
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In Q!i, a2, as and (3 are 4x4 complex matrices, whose standard form (in 
2x2 blocks) is 

where I = ( r, i ) ^'^d Ufc are the Pauli matrices: 



1\ /O -A 1 



1 



The magnetic Dirac operator without the Coulomb potential has essential 
spectrum (— oo, —1] U [1, oo) and no eigenvalues in the gap (—1, 1). For v e (0, 1) 
the Hamiltonian has the same essential spectrum and eigenvalues in the gap. 
The ground state energy \i{v,B) is defined as the smallest eigenvalue in the 
gap. 

As the field gets large enough, one expects that the ground state energy of the 
Dirac operator decreases and eventually penetrates the lower continuum. The 
implication of this for a second quantiz ed model is that electron-positron pair 
creation comes into the picture |18II20| . The intuition for that can be gleaned 
from the Pauli equation, where the magnetic field tends to lower the energy 
because of the spin. It is therefore reasonable to define the critical field strength 
B{v) as the supremum of the positive B's for which Ai(i^, b) is in the gap (—1, 1) 
for all b S (0, -B). As a function of v, \i{v,B) is non-increasing. Hence the 
function B{i') is also non-increasing. 

One of our goals is to give estimates on this critical field as a function of 
the nuclear charge. Our first result, proved in Section 2 is that this critical field 
exists and we give some rough estimates in terms of v: For some C > 0, 

As a corollary we get the noteworthy result that as v ^ 1 the critical field B{v) 
stays strictly positive. This is somewhat remarkable, since in the case without 
magnetic field the ground state energy as a function of v tends to as — > 1 
but with an infinite slope. Thus, one might expect very large variations of the 
eigenvalue at ;y = 1 as the magnetic field is turned on, in particular one might 
naively expect that the ground state energy leaves the gap for small fields B. 
This is not the case. Moreover, since the hydrogenic Hamiltonian ceases to be 
selfadjoint at = 1 it is hard to visualize how one might arrive at such estimates 
using standard perturbation theory. 

Section 3 is devoted to the asymptotics of B{v) as ^ 0. We define the 

notion of lowest relativistic Landau level which leads to a one dirnensional ef- 
fective theory. This effective theory can be analyzed in great detail and allows 

to calculate the ground state energy X^iv^B) of the magnetic Dirac-Coulomb 

equation (^) in the lowest relativistic Landau level. It is given by the variational 

problem 

\'~{v,B):^ inf A^[/,i/,B], 



where A = A''' [/, v, B] is defined by 

a/ |/(z)pdz= / {—l^^^^+[l-j,a^{z))\f[z)\^]dz 



and 



+00 gg-iS, 
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The point here is that for B not too sniaU and i' not too large (the precise 
bounds are given in Theorem 0)), 

The one dimensional Af(i^, -B) problem, although not trivial, is simpler to 
calculate than the Xi{v, B) problem. As a result, in the limit as — s- 0, this new 
theory yields the first term in the asymptotics of the logarithm of the critical 
field. In particular we have the following result, 

lim i'\og{B{iy)) = TT . 



From a methodological point of view, the ground state energy of the Dirac 
operator is not given by a minimum problem for the corresponding Rayleigh quo- 
tient, but it is a min-max in the sense that one decomposes the whole Hilbert 
space Ti. = TLi © TL2 , maximizes the energy over functions in I-L2 and then mini- 
mizes over non-zero functions in Tii. While the choice of these Hilbert spaces is 
not arbitrary, there is some flexibility in choosing them, see IT, 10 ,5 . For certain 
choices, the maximization problem can be worked out almost explicitly leading 
to a new energy functional for which the ground state energy Ai(i/, i?) is the 
minimum. In this sense, the ground state energy of the Dirac operator appears 
as a minimum of a well defined functional. Both variational characterizations, 
the min-max and the min, are of course equivalent, and our approach depends 
on the interplay between the two. 

Our results are different from the work of [2] which considered the non- 
relativistic hydrogen atom and worked out the asymptotics of the ground state 
energy as B ^ 00 for every v > 0. In our case, however, v has to stay in the 
interval [0, 1) in order that the operator can be defined as a selfadjoint operator. 
Further, the critical field is always finite and we are interested in estimating it 
as a function of v. The similarity with j2| comes as we let v ^ 0, since then the 
critical field tends to infinity but the estimates are not the same. 

While the mathematical methods are the main point of this paper, let us 
make a few additional remarks about its physical motivation. Spontaneous pair 
creation in strong external fields, although never experimentally confirmed, has 
been analyzed by Nenciu [T7I[TH| . In ^7] it was conjectured that by adiabatically 
switching the potential on and off, there is spontaneous pair creation provided 
some of the eigenvalues emerging from the negative spectrum crossed eigenvalues 
emerging from the positive spectrum. This conjecture was partly proved in 
and 12 Oj . Since such a crossover occurs in the Dirac hydrogenic atom with a strong 
magnetic field, it is natural to try to estimate the strength of the magnetic fields 
for which this crossing phenomenon occurs. 

Note that the unit in which we measure the magnetic field is huge, about 
4.4 X 10^ Tesla. Sources of gigantic magnetic fields are neutron stars that can 
carry magnetic fields of about 10^ Tesla. Fields of 10^^ Tesla for a neutron star in 
its gestation and in magnetars are expected, and there is speculation that fields 
of up to 10^2 Tesla may exist in the interior of a magnetar. There is a considerable 
literature in this area and an entertaining introduction can be found in \Ti\ . 

Further, it is expected that near the surface of a neutron star atoms persist up 
to about Z — AO. We show that the critical field at Z = 40 must be larger than 
4.1 X 10^" Tesla, and preliminary calculations using numerical methods based 
on Landau levels yield an estimated value of about 2.5 x 10^^ Tesla. Although 
improvements on these estimates are currently under investigation we believe it 
is unlikely that they will yield relevant values for the magnetic field strength. 
For elements with higher Z, the values for the critical field are much lower. In 
the case of Uranium {Z = 92), they are sandwiched between 7.8 x 10^ Tesla and 
an estimated value (using Landau levels) of 4.6 x 10^^ Tesla. 

Speculations that large magnetic fields facilitate the creation of electron - 
positron pairs are not new in the physics and astrophysics literature. Clearly, 
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the Dirac operator coupled to a magnetic field but without electrostatic potential 
has a gap of 2rnc^ independent of the magnetic field. It was pointed out in 
I19| that the anomalous magnetic moment narrows the gap, i.e., it decreases 
the energy needed for pair production. In lowest nontrivial order the anomalous 
magnetic energy is proportional to the magnetic field which leads indeed to a 
narrowing of the gap; in fact the gap closes at a field strength of about 4 x 10^^ 
Tesla. It was observed in [T^, however, that the anomalous magnetic energy 
depends in a non linear fashion on the external field. Further it is shown that 
even at field strengths of 10^^ Tesla the gap narrows only a tiny bit, irrelevant 
for pair production. For a review of these issues the reader may consult 8 . Our 
contribution is to take into account simultaneously the magnetic field and the 
Coulomb singularity, in which case no explicit or simple calculations are possible. 

Of course our analysis only deals with a single electron and a fixed nucleus. 
A description of the non-relativistic many electron atom under such extreme 
situations has been given in 9,15,16]. The authors study various limits as the 
nuclear charge and the magnetic field strength gets large and determine the 
shape of the atom in these limits. In non-relativistic physics the natural scale 

for the magnetic field is w 2.4 x 10^ Tesla, much smaller than the ones 

under considerations in our paper. 

As we have mentioned before, for small Z the critical magnetic field is of 
the order of e^^ and hence non relativistic physics is sufficient to explain the 
rough shapes and sizes of atoms even at very high field strengths. It may very 
well be, however, that for heavy elements and very large fields qualitatively new 
effects appear that cannot be understood on the basis of non relativistic physics 
alone. Should such effects occur, then it could make sense to treat the many 
body relativistic electron problem using the Dirac - Fock approximation. 



2. Ground state and critical magnetic field 

In this section, we set some notations, establish basic properties and prove Es- 
timate (j21l on the critical magnetic field. 

2.1. Min-max characterization of the ground state energy. The eigenvalue equa- 
tion for the Hamiltonian 

Hb^-^^ = X4' (3) 
\x\ 

is an equation for four complex functions. It is convenient to split ^ as 




where (f),x ^ L'^{M.^; C^) are the upper and lower components. Written in terms 
of (j) and Xj ® is given by 

^sX + 0-/^0 = A0 , (4) 
\x\ 

PB(t> - X - 1^ X = Ax . (5) 

Here Pb denotes the operator 

Pb ■.= -icr- {V -iAeix)) , 

where 



Ab{x) 
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is the magnetic potential associated with the constant magnetic field 




Using |(SJ) we can eliminate the lower component % in Taking then the 
inner product with cj) we get 



J[(j), A, I/, = , (6) 



where 



JK3 \ 



Thus, we see that that for any eigenvalue A G (— 1, 1) of the corresponding 
eigenvector leads to a solution of lO- 

Reciprocally, the functional J can be used to characterize the eigenvalues. For 
this purpose, a few definitions are useful. The functional J[0, A, i^, B] is defined 
for any B &m+,v & (0, 1), A > -1 and </) e C^(R3^C2). Further, in order that 
(jSJ makes sense we introduce the set 

A{v,B) :={<^e Co°°(M3) : UU^K^) = 1, 

A ^ J[0, A, V, B] changes sign in (—1, +oo)} . 

Note that this set might be a priori empty. Finally, since the function J is 
decreasing in A, we define A — X[(j), v, B] to be either the unique solution to 

J[(j3,\v,B] =Oif0e A{v,B) , 
or \[(f), V, B] = -1 if J[(f), -1, V, B] < 0. 
Theorem 1. Let B eR+ and v e (0, 1). // 

-1 < inf X\(t), i',B] < 1 , 

this infimum is achieved and 

Xi{v,B):^ inf A[(/>,j/,B] 

<t>eA(v,B) 

is the lowest eigenvalue of Hg — j^ in the gap of its continuous spectrum, (—1,1). 

Proof. This proposition is a consequence of Theorem 3.1 in 0. The essential 
assumptions of this theorem arc: 

i) The selfadjointness of Hb — v \ ■ which is proved in the appendix. It is 
crucial here that < < 1. 

ii) The existence of a direct decomposition of L^(M'^;C^) as the sum of two 
subspaces Hi © H-i such that 

(x,{HB-y\-\-^)x) 

a-i := sup 

. f + [HB-y\-\-^){x + y)) 
< ci :— mt sup . 
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Set b ini (Tcss{Hb ^ \ ■ \ ^) H (02, +00). li ci < b then ci is the lowest 
eigenvalue of Hb — v \ ■ in the interval (02, h). In the present case we choose 
the decomposition 

, (<i>\ f(t,\ /O^ 



based on the upper and lower components of the four components spinor ip. It 
is easy to see that a-z = —1. Furthermore the essential spectrum of Hb — | • 
is (— c», —1] U [1, +00) independently of B, see Hence 6 = 1. It remains to 
calculate the supremum in the definition of ci as a function of x — (p). Note 
that the Rayleigh quotient in the definition of ci is strictly concave in y — (°) . 
Therefore the supremum is uniquely achieved by 



and its value is A[0, ly, B], that is, 



A[0, ly, B] = sup 

xeCo-(R3,C2),^=(*) 



(V',V') 



□ 



Remarks: 1) The eigenvalue Xi{i', B) can be characterized either as the min- 
imum of the functional X[(j), B\ or as a min-max level of Hb — v \ ■ . Both 
characterizations will be useful in the sequel of this paper. 

2) Under the assumptions of Theorem^ we have 

J[(l),\,iy,B]>0 V(/)e C5"(]R^C2) 

for any A < Xi{i',B). The eigenvalue Xi{h',B) can therefore be interpreted as 
the best constant in the above inequality. 

3) When Ai (1^,-8) is equal to —1, it belongs to the continuous spectrum and it 
is not necessarily an eigenvalue of Hb — v \ ■ . 



2.2. Basic properties of the ground state energy. 

Proposition 1. For all B > 0, the function v ^ Xi{h',B) is monotone nonin- 
creasing on (0,1). 

The proof is left to the reader. It is a consequence of the definition of J[(/), A, ly, B]. 

Proposition 2. For all B > 0, the function v — > X\[v^B) is continuous in the 
interval (0,1) as long as Xi{i', B) G (—1,1). 

Proof. By Theorem^ if Ai(i^, B) e (—1, 1) there exists a function (j)^, such that 
J[(f)^, Xii^v, B), v, B] — 0. For any sequence converging to ly, the upper 

semi-continuity of — > Ai(i^, B) holds: 

limsup Ai(f„, B) < limsup A[0i/, f„, -B] = Ai(j^, B) . 

n — *-t-oo n — >+oo 

If Vn < V, then A„ := Ai(t'„,i?) > Xi{v,B) and {A„}„ converges to A. Consider 
therefore a {t'n}„ converging to v from above. We have to face two cases: 
First case: Xn > —1 for all n e N. Since J[4)u,Xi{v,B),Vn,B] < 0, we know 
that A„ < Xi{i',B). Consider the corresponding eigenfunctions ipm such that 
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J[(t>m An, i^m B] ~ 0, where denotes the upper component of ^„ and assume 
that ||0n||L2(R3) = 1. By Theorem ^ we have 

l^-S*^"!^ , /I \ Wa. \i\ _ /" ^"1- |2 .3, 



1 + A„ + 



— + (1 - A„)|0„|^ d-'x = / ^ |</.„|^ d'x 
R / 



Assume that A := hminf„^_|_oo -B) < Xi{h',B). Up to the extraction of 

a subsequence, assume further that {Ai(i^„, i?)}„gN converges to some value in 
{—l,Xi{i^,B)) and choose A G (yl, Ai(^, i?)). For n large enough, Ai(^'„,i3) < A 
and 

- + (1 - mnA d^x< I M ■ (7) 



Second case: \i{v\B) = —1 for all > v. We choose A G (— 1, Ai(iy, B)) and 
find a {(^„}„ such that ||(/)n||L2(R3) = 1 and J[0„, A, iy„, S] < for n large: iQ 
also holds. 

Using the monotonicity of the {l'n\n^ which implies the monotonicity of the 
{A„}„ by Proposition n and the fact that in both cases, v G (0,1) and A G 
(— 1, 1), we deduce from |(7J) a uniform bound for the functions 

sup / (\x\ |Pb0„P + ^) d^'x < +O0 . (8) 

The proof goes as follows. It is sufficient to prove that Jj^a |a;|~^ |0„pd'^a:: is 
uniformly bounded. Let x be a smooth truncation function such that x(r) = 1 
if r G [0, 1), x{r) = if r > 2, and < x < 1. Since 

with (j)n = x{\R^^ • I) 't'n, it is therefore sufficient to prove that J^s \x\^^ l^nP d'^x 
is uniformly bounded, for some i? > 0, eventually small. Using the estimate 

^2 > + _ 

- 1+e e ' 
we get the following lower bound 



. [ \PbK? ,3 C ~ 2 

-L (l + s)(l + A+ft)'""""^"""^"^ 

where C is a constant which depends on ||x'lli=o(K+)j B and R. Next, with the 
same type of arguments, we can write 



1 + e £ 



> 



1 + e 



8 



Jean Dolbeault, Maria J. Esteban, Michael Loss 



Collecting these estimates, this gives 

■£x<C{e,R,x)+ I ^l^n^'d^'x 



(l + e)2(l + A + |^) - ^ ' ^3 \x 

Because 0„ has a compact support in the ball of radius 2R, if 6 = 6{R) > 
(1 + X)R/vn at least for n large enough, then 

1 \x\ 

= > , ' VxeS(0,i?) 

1 + A + 1^ - J^n (1 + <5) 



so that 

1 



JR3 JR3 FI 



On the other hand, according to 0111, 

\x\ \a-V^n\''d'x > I ^l^npd^rr . 

Jr3 \x\ 

This provides a uniform upper bound on Jjja |a;|~^ d^x if £ and 5 are chosen 
small enough in order that 

1 2 

> K 



(l + £)2(l + ,5) 



for n large. This can always be done since Vn converges to £ (0, 1) and 5{R) 
can be taken as small as desired for i? > sufficiently small. This concludes the 
proof of ©. 



Summarizing, we obtain that 



J 



<t>n,\{\ + \i{v.B)),v,B 



<o , 



for n large enough: hence A [(/)„, t^, S] < ^{\ + \i{v,B)) < Xi{i',B), a contra- 
diction. □ 

Consider now the effect of a scaling on J. 

Lemma 1. Let B >0, X> -1, 9 > and <j>g{x) := e^/'^(p{ex) for any x e M^. 
Then 

yA,2j9ix)^0'^/^[V(j}iex)-iABiex)(i}iex)] , 

and for any a G M, e (0, 1), 



j[ci>e,x,0'^^,e^B]= I {e \ l^^-^l'^^ +(l-A)|^p-g^H^) d'x 



1 + A- 



-R- 



(9) 



Using this scaling, we prove some properties enjoyed by the function \i{v,B). 
Take A = X[4>e, v, O'^B] > -1, 6* > 1 and a = in 10. With 1 + A = 6'(1 + fi), 

= J[4>e,X,iy, e^B] 



1 + Ai 



\x\ 



i0n d^x-f 2(1-0) / 

/ JR3 
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Assuming that ||^!>||l2(r3) = 1, we get 

J[</),M,i^,S] =2^ >0 

and thus, 

A ^ — 1 

\[ci),v,B\> n=--—^ . 

On the other hand, -^J[(j), fi,v,B] < —1, so that 

^ Q ^ 

X[ct>,i^,B]<^i+J[ct),n,iy,B] = - + -^ . 
Summarizing, we have the estimate 

6 9 6 9 

The above estimate, which holds provided X[(j)0, u,9'^B] > —1 is equivalent to 

*i^m_^<,(,,,.,„,Bi<*f!5! + ^ v«>i (10) 

under the condition \[(f), v, 6'^B] > —1. As a consequence, we have the following 
result. 

Proposition 3. For all v £ (0, 1), the junction B ^ \\(y, B) is continuous as 
long as it takes its values in (— l,+oo). Moreover 

M'^,0'B) e-1 X,{u,6^B) 6-1 

— 9 ^'^"'^^ - — 6 — + -r ' ^^^^ 

if\i{v,B) e (— l,+oo) and 6 e T^xi^irB)) • ^ consequence, B i— > \\{v,B) 
is Lipschitz continuous for any u e (0,1) and B > such that Xi{i',B) e 
(-l,+oo); 

Ai - 1 ^ aAi ^ Ai + 1 
2B - dB - 2B ' 

Proof. Choose a e {-l,Xi{u,B)) and take any (f) G C^{M.^,C'^)- Since 

^Jlcl,,X,„,B]<-l , 
an integration on the interval [a, X[<j), v, B]] shows that 

r -1 A=A[0,i/,B] 

-J[(j),a,i',B]= J[(j),X,i',B] < -X[(j),i',B]+a 

L J A=a 

where the first equality holds by definition of A[^, B] , i.e. J[(j), A[^, u, B] ,p,B] = 
0. As a consequence, 

J[(/>, a, u, B] > X[(l), I/, B] - a > . 

The function 6 i— > J[<p, a, v, 6^B] is continuous, so only two cases are possible: 
First case: 

J[(t), a, V, 9^B] > V 6» > 1 , 
Second case: there exists a constant 6 = 6{a, (p) > 1 such that 

(i) J[0, a, V, e^B] > for any 6 G (1, 9), 

(ii) J[(j), a, v, 6^B] = or, equivalently, X[4>, v, 6^B] = a. 
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In the second case, by (i), we know that X[(j>, ly, 9^B] > a > — 1 for any 6 E (1, 6) 
and so ifTT))! apphes: 

In the hmit case 6 — 6, we get 

e\i{v,B) + 1-9 < X[<j3,v,pB] = a , 
which gives the estimate 

0> , =:e*(a) . 

- 1 - Xi{iy,B) ^ ' 



Thus the inequahty 



Xi{v,B) < 9X[(Pi/g,iy,B] < X[<p,jy,0^B]+9-l 



holds for any 9 e [0, 9*{a)] and for any </> g C^(M^, C^), which proves the r.h.s. 
inequahty in ((TT|l by letting a ^ — 1: 

2 

lini 9* (a) = 



-1 ^ ' l-Xiiiy,B) ■ 

The l.h.s. inequality is obtained in the same manner. □ 

With the appropriate test functions one can prove that Ai {v, B) is always 
below — 1 for i? large. We recall that Ai(i^, B) — —1 means that if J[(j), ~l,iy,B] < 
for any 4> G B). Let us give some details. 

Proposition 4. Let v £ (0, 1). Then for B large enough, Ai(zy, B) < and there 
exists B* > such that Xi{iy, B) = —1 for any B > B* . 

Proof. Let us consider B > and the trial function 



^0 
where 

^ y 2tt \ J 

and / e Cg°(M,R) is such that / = 1 for |x| < S, S small but fixed, and 
||/||j^2(B) = 1- Note that G Ker{PB + icr^ dx^) and so, 

Pb0 ^ -i (738x3 (/> , 
Moreover, the state cj) is normalized in i^(R'^). With r = \x\, we can define 



G 



B\ 



(- IPb^/'P - - d^x (12) 
\i/ r J 



and compute 
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Using ^/B-Ks^ + ja^gP < B-1/2 g + 1^:31 and 



// 



> log(52i3) 
4ve 



for B > 1, we can therefore bound Gb[0] from above by 

Gs[0] < ^ + C2i^ -CsiylogB . 

where Cj, i — 1, 2, 3, are positive constants which depend only on /. For 
B > 1 large enough, Gsl^] + 2||(/)|p < and Ai(z^, S) = —1, since in this 
case J[(/),- 1,1/, B] < 0. □ 



2.3. The critical magnetic field. Proposition^niotivates the following definition. 
Definition 1. Let v G (0, 1). Wc define the critical magnetic field as 

B{v) := inf |b > : hm Xi(y, b) = -l| . 

Corollary 1. For all v £ (0, 1), \i{v,B) < 1 for any B e (0,B(i^)) . 

Proof. For S = we have Ai(i/, 0) = \/\ — v"^ < 1. Given B > 0, small, 
by continuity of £? i-^ Xi{i',B) we know that Xi{i/,B) S (0,1). Let us con- 
sider e (1, ^jB{v)/B) such that -1 < \i{v,e'^B) < 0. This is made possible 
by Propositions 121 and ^ Then, by Proposition O 

Ai(z^,B) < ^ < 1 . 

□ 

The computations of Proposition 0] show the existence of a constant C3 > 

such that B{i') < e^^^" , for all i> £ (0,1). This estimate can be made more 
precise for any v not too small: 

Theorem 2. For all v G (0, 1), there exists a constant C > such that 

The proof of this theorem uses Proposition^] Otherwise it is splitted in two 
partial results stated in Propositions |5l and |2| 

Notice that there is a big gap between these lower and upper estimates when 
1^ is small. To try to better understand this problem, in the next section, we will 
analyze the case when B is large, proving that the 3d definition of B{i') is actu- 
ally asymptotically equivalent to a Id problem related to the lowest relativistic 
Landau level. More precisely we will prove that when 1/ is small, B is not too 
small and Xi{i', B) G (—1,1), the eigenvalue associated with Xi{i',B) G (—1,1) 
is almost equal to the corresponding eigenvalue in the lowest relativistic Lan- 
dau level class of functions, see Theorem 0] We will also establish that B^i/) 
behaves in the limit v like the upper bound in Theorem |21 and obtain the 
corresponding value of C, see Thcorcm|Sl 

Our first partial result is the following 
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Proposition 5. For any v e (yS/S, 1), \J B{v) < ^ 



3 



Proof. Consider the trial function t/j — (g) where 



B 



3/4 



is like the one chosen in the proof of Proposition^ with /{x^) — {^)^^'^ e ^'^l/'^. 
Here, with the notation r := |a;|, we find 



= (27r)"5 B 



— / e-'"'' '"^ r^ dr cos^ sin 9 dO - A n 



''"/^rdr 



3 1 



TT 

4 TT 1/ 



^/'^r^dr - Am 



1'^ r dr 



If 1/2 e (2/3,1) and ^/B > ^^0^, then Gi3[0] < -2 = -||0||2 and so Ai (i/, S) = 
— 1, which proves the Proposition. □ 

Proposition0]shows that for > ■\/2/3 and i? large, B) possibly ceases 
to be an eigenvalue of the operator Hb — v\ -\~^. This can be interpreted by saying 
that for strong magnetic fields, the Coulomb potential does not stabilize the 
electron. At some level, electron-positron pairs could appear and then Quantum 
Field Theory (or QED) becomes unavoidable for a correct description of the 
electron dynamics, see PH] . 

Proposition 6. For given v E (0, 1) and B > Q, the function 9 i-^ Xi{9^^h', 9^B) 
is monotone nondecreasing as long as it takes its values in (—1, 1) and i>/9 G 
(0,1). 

Proof Take a = -1 in ®: 

J [<^e, A, 0- V, 9'B] = {9' - 1) / + J[<p, A, V, B] 

Jr3 1 + a + 

so that for < 1, 

J [(j>e,X, 9-^v, 9^B] < J[4>, A, i^, B] 

at least for 1 — 9 > 0, small, so that (f>g e A(9^^u,9^B) for (f> = (f>,j such that 
A[(/i^, V, B] — Ai(i/, B). This proves that 



\i{9-^i',9^B) < Xi{i^,B) (13) 

for 1 — 9 > 0, small. By continuation, the property holds as long as the assump- 
tions of Proposition are satisfied. The case 9 > 1 follows by multiplying 9 v 
and 9~^B by respectively 9~^ and 6'^. □ 

Corollary 2. There exists a positive constant A such that 



B{v)> 



A 



as v \ 
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Proof. Let {lyo.Bo) be such that B{vq) > Bq, i.e., Ai(j'o,-Bo) > ~1 and take 
ly e (0, i^o), e ^ v/vq e (0, 1), B ^ e-^Bo in 



-1 < Xi{iyo,Bo) < \i{v,B) = Xii^Bo/Bi^o.B) 
By Proposition ^ this inequahty can be extended to 



Xi{^Bo/Biyo,B) < \i{v,B) V G (0, ^B^/Bvo) 
This amounts to say that 

< < ^JB^/B =^ B{v) > Bo 4 , 
which proves the result with A = Bq i/q . 



□ 



The constant A can be made more precise. The remainder of this section is 
devoted to the following improvement of Corollary |21 



Proposition 7. For all v ^ (0, 1), 



V 



m? V Jr3 \x\ 

In particular this implies that B{v) > 
Proof. Scaling the function <f) according to 

cj^B ■.= b^'^4>{b''^x) 

preserves the norm, and yields 

Gb[4>b]^^Gi[4>] , 

where Gb has been defined in p2fl . Obviously it is sufficient to find a good 
estimate on the functional Gi. 

Let us collect some preliminary observations. Recalling that the angular mo- 
mentum vector L is given by 

L = -iV A X , 



a simple calculation shows that 
tJ- V = 

with r = Ixl and 9^ = - • V. We also recall that 



1 



-i Ai(x) ■ a — — (cr ■ — ] (a ■ q{x)) where q(x) — — 
\ r J 2r 



-X3 Xi 
-Xz X2 



and 



i Pi = a ■ \7 ~ i Ai{x) ■ a 
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SO that we can expand as 



dr a ■ L — a ■ q{x) 



\drct>\'' + ^\u-Lcj>\'' + \q 



2 I J,|2 



dr{(f), a ■ q(f)) + {(j),a ■ (drq) (f)) 



dr{4>, a ■ L (j)) + {a • L (j), a ■ q (p) + {a ■ q (j), a ■ L I 



As a last preliminary remark, we notice that r dr-q = q. 

Since the vector potential grows linearly, we localize the problem near the 
origin. To this end consider the function 



t{r) 



1 ]ir<R, 
R/riir>R . 
Since t{r) < 1 and f < 1 we get the lower bound 

t{r) r 



Gi[4>] > 



> 



|2 j3 



y( f t(r)r\PiM^d^x- f - IM^ d^x^ = ( K[S - [ -IS'^d'^x 



where the kinetic part is defined by K[(j)] := t{r) r \Pi(j)\^ d^x and satisfies 

K[cf>] 



= / t{r)r 
+ I t{r) 

JR3 



\drci>\' + ^\a-Lcl>\' + \q\^\ct>\^ 



d'^x 



- r dr{(t>,(J-q4>) + r {<f),a-{drq) (fy - dr{(j),a-L(j)) 



d^x 



+ / t{r) {a-L(j),a-q(j)) + {a-q<j),a-L( 
An integration by parts in the r variable yields 

m 



d^x . 



t{r) r 
+ / t{r) 

JR3 



d^x 



4:(^,a-q^) + - {(j},(j-L<f)) + {a-L(f),a-q^) + {a-q^,a-L(t)) 
r 



d^x 



+ / t\r) r{(j),a-q^) + {(j>,a-L> 



R3 



d-'x 



where we have also used that rdrq = q. Consider now the region of integration 
where r < R and denote the corresponding expression by -f^i There the 
derivative of t{r) vanishes and hence collecting terms we find 



/ tr 

Jr<R 



-{a-L+l) + a-q 
r 



+ 2 



d'^x 



r<R 



,a ■ qtj)} d^x . 
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At this point wc have decoupled the derivatives with respect to r from the 
magnetic field or, to be precise, from q. The problem is that the angular mo- 
mentum is still coupled to the magnetic field. Obviously 



(O- • L + 1) + £7 • 5 



> 



\a-q( 



Further, since the eigenvalues of a ■ L + 1 are given by ±1, ±2 . . . , 

II ((T-L-|-l)(^ 11^2(5-2) > ||0||i2(S2) , 

and we have that 



-{a-L + l) + cr-q 
r 



1 



L2(52) 





> 


1 r- 










r ~2. 



since \q{r)\ < r/2. For r < R < -\/2, the factor [i — 5] is nonnegative. Since 
ir ~ i]^ ~ 7^ = T ~ since 2 \q{r)\ < r, we obtain the lower bound 



Ki[<P] > 



r<R 



-2 



The function 2^ — 2r is a decreasing function on the interval [O, v^] . Hence 

(fx , i{R<V2 . 



Ki[(t)] > [ r\dr(t)fd^x- 

Jr<R 



i?3 



r<R 



Next, we look at the contribution to K\(jj\ of the region where r > R, which 
we denote by -fsr2 [</>], 



= t / t{r)r 

Jr>R 



d^x 



+ 



4:{(j),a-q(j)) + - (4>,a-L(j)) + {a-L(j),a-q(j)) + (a-q(j),a-L< 
r 



I t{r) 

Jr>R 

[ \ r{cl>,a-qct>) + {<!>, a- L<f>) 



using the fact that t' = —t/r. Collecting the terms, we get 

1 , , ,,0 , .0, ,,ol^3^ 



= / t{r)r 

Jr>R 

+ I t{r) 

Jr>R 



i5.</.i^+-ic7-i0|2+i9ri</)i^ 



1 



S((j),a-q(l)} + - ((j),a-L4>) + {a-L(j),a-q(j)} + (a-q(j),a-L( 



This can be rewritten as 



\<t>]> I t{r)', 

Jr>R 



a.L+-_ 



+ cr ■ q 



4^2 



d-'x . 



d'^x 



+ 2 t{r){(j),a-q(j))d:^x . 

Jr>R 
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Finally we get 

K2[4>] > I t{r) r |5,<^p d?x-^ ( |<^|2 d?x -2R [ - \q\ |<^|2 d'x , 

Jr>R Jr>R Jr>R ^ 

and, using \q\/r < 1/2, 



K2[(l>] > / t{r)r\dr<l>\^(Px- 

Jr>R 



R + 



AR 



\(^Yd^X . 



r>R 



Thus we can estimate K[(f)] = Ki[<p] + K2[(f)] as follows: 



K[(f)\ > [ t{r)r\dr(f>fd^x+ 
Observe that generally 



i?3 



-2R 



r<R 



R+ 



AR 



\cl)\U-'x 



r>R 



dr(t> H I 

r 



< / t{r) r 

= [ t{r) r \dr(i)\^ d^x - [ t{r) - d^x - [ t'{r) \<f>f S 



X . 



Since t'{r) = and t{r) = 1 on [0,R), and t'{r) = -t{r)/r on (i?,oo), the 
following estimate holds 

/ t{r)r\dr^\U^x- [ -\<l>\^d^x>- [ ^^^\cl>\U^x+ [ t' {r) \ct>\U^ x 

= -/ hcl>\'d'x>-U \ct>\'d'x 

Jr>R ' ^ Jr>R 



and hence 



m- [ -\^\'d' 



X > 



R-^ 

--2R 



Jr<R 



J2 j3 



d-'x . 



r>R 



Optimizing on i? e (0, V^], i.e. using 



max mm 
fle(o,%/2) 



i?3 



2R 



R 



AR 



R 



AR 



\R=V5/2 



= -V5 



we get 



Gi [(/.] > u (K[(t>] - I - d^x] > -z/ \/5 / d^x . 
Hence the condition G'_b[(/'] = •\/BGi[0] > —2 entails that 



B{y)> 



□ 
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3. Asymptotics for the critical magnetic field 

In the large magnetic field limit, the upper component of the eigenfunction cor- 
responding to the lowest energy levels in the gap of Dirac operator with magnetic 
field Hb — I • is expected to behave like the eigenfunctions associated to the 
lowest levels of the Landau operator 

Lb ■■=-i CTi - i<^2dx2~cr ■ Ab{x) , 

which can also we written as Lb ^ Pb+ i (^z dxa ot Lb ^ — i (dxi + i BX2/2) ui — 
i {dx2''iBxi/2) (72. The goal of this section is to compare the lowest energy levels 
of Hb — v \ ■ with its lowest energy levels on a space generated by the lowest 
energy levels of Lb- The asymptotic analysis for the small coupling limit v ^ 0^ 
is not that simple because the Landau levels are not stable under the action of the 
kinetic part of the Dirac Hamiltonian. The way out is to choose a representation 
of Hb — v \ ■ that diagonalizes the kinetic energy in the Dirac Hamiltonian 
and to project both upper and lower components on the lowest Landau levels. 



3.1. Projection on Landau levels. To start with, we observe that 

p2 _ r 2 q2 
^B — ^B '~'x3 

and summarize the basic properties of the lowest energy levels of L^. 

Lemma 2. Section 7.1.3] The operator Lb in L^(R^, C'^) has discrete spec- 
trum {2nB : n £ N}, each eigenvalue being infinitely degenerate. Moreover the 
kernel of this operator, that is, the eigenspace corresponding to the eigenvalue 0, 
is the set generated by the L^ -normalized functions 

r>(e+l)/2 /n\ 

^'■=- .{x2+ix,Ye-''' ], £eN, s'^xl-^xj. 



^/¥K2n\ vi 

Next we diagonalize the free magnetic Dirac Hamiltonian. First we write it 
in the form 

Pb\_ 177^ (R Q 



where R and Q are operators acting on 2 spinors, given by 



and satisfy the relation 

i?2 + Q2 ^ I . 

The matrix 

'R Q 

Q -R 

is a reflection matrix and hence has eigenvalues 1 and —1. It can be diagonalized 
using the matrix 

1 ( Q 
^20— R) \}-R Q 

The operator defined by U is unitary and such that 



U*KbU= ( ^^ + ^B 



-v/TTp] 
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The potential V = is transformed into the nonnegative operator 

P:=U*VU= (^^, ^ 

Here and from now on, we will omit I whenever it is multiplied by a scalar valued 
function. If wc denote by Z any 4-spinor and decompose it as 



Z 



where X and Y are the upper and lower components, in the new representation, 
the full magnetic Dirac Hamiltonian takes the form 



The Dirac energy for an electronic wave function Z in the electromagnetic 
potential (V, A) is now 

£,[Z] := IC[Z]-i^{Z,PZ) , 
IC[Z] := (Z, U*Kb UZ) = (X, Vl + Pix) - (y, ^/TTp^Y^ . 

As we shall see below, in the new representation, restricting the upper and lower 
components X and Y to the lowest Landau levels makes sense for studying the 
regime of asymptotically large B. The price we pay for that is that all quantities 
like R, Q, U, P... depend on B. Denote by Be the projector on the lowest 
Landau level, whose image is generated by all functions 

and define 77^ ;= I — He- Notice that He commutes with Lb- With the above 
notations, for aU ^ e L'^{m.^,€?), 

we have that 



+ ^ j = (^^^^' + nci ) + [nii., VI + Pb m ) • (i4) 

Next, we decompose any Z E (iy^(M^, C))"* as 

Z = B Z + B" Z , 

where 

■ ' Be) ' V ^ ^£ 



3.2. Main estimates. From 1141) . it follows that 

K.[Z] =K.[B Z]+K.[B'' Z] . 
Since the operator P is nonnegative, we also have 

{Z,PZ) <{l + ^) {BcZ,PBZ)+ + {B^Z,PB'Z) 

{Z,PZ)> {BcZ,PBZ)+ (l-^) {B' Z,PB''Z) 
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This simply follows from the identities 
(a + 6)2 < + + 2\ab\^ inf 
(a + bf > + - 2\ab\ = sup 

u>0 



The above remarks prove the next proposition. 
Proposition 8. For all Z e Cg= {M.^ ,C'^], 

[n^z] < £,[z] < £,_,.,^[nz] + £„_^[wz] . 

The following result will allow us to get rid of the higher Landau levels when 
looking for the ground state energy, i.e. of the term £yj^^ [n"^ Z] in Proposi- 
tionlHI Consider D e (0, 1) such that 

2{D + V^) = 2-V2, 

i.e. V « 0.056, and for any v E (0, D), define 

d{6) -.^ {l~2d)V2-26 , d±{i^) := did±{v)) with S±{v) := ± ly . 

We have 

d{S) > S <1- V2/2 , 

d±{iy) > if 1^ < D . 



Proposition 9. Let B > and S £ (1 - \/2/2). For any Z 

X,Y e L^{R^,C^) 



£s [n^z] > d{S) VB \\n^x 



£-5 [n" Z] < -d{S) VB II 77^ Y 
Proof. An elementary computation shows that 

P| = {V -iAefl + a-B 
Using the diamagnetic inequality (see I ), 



Hardy's inequality, 



viv^l 



dx > 



1 



V|^| 



dx 



dx > — - — dx , 



and the nonnegativity of ct • B + _B I, we get 

+ (1 + B) I = (V - iAef I + o- • B + (1 + B) > 
Since the square root is operator monotone, we have 



1 1 

4 bp 



1 1 



VB + Ji + P|> JP| + (1 + 5)1 > J- -^ + 1 



4 x 
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Now, for any 5 > 0, 



> (1 - 26) ^yi + P'i -2SVB . 



On the range of 77^ the operator + P^> ^/T+L^ is bounded from below 
by VI + 2B. Hence 



/l + P| - — > (1 - 2,5) V1 + 2B - 2,5 /B , 
which is equivalent to 



/l + P|-5P> (l-2,5)Vl + 2P-2J\/P . 
Since 

(1 - 2,5) VTT2B -2SVb ^ d{5) \/P 

as _B — > CX3, and 

inf [(1 - 2S) VTT2B -2S^/B]>0 

B>0 

if d{S) > 0, the right hand side is positive for any field strength provided that 
0<S <1- V2/2. □ 



3.3. The restricted problem. Next we prove that the ground state energy Ai (v, B) 
is comparable with the one obtained by restricting it to the lowest Landau level, 
both in the upper and in the lower components of the wave function, provided 
that the Coulomb potential is slightly modified. By a result similar to Theorem 
□ (also see Theorem 3.1 in in case B = 0), for all < P < B{v), v e (0, 1), 
Ai(i/, B) is characterized as 

\i{v,B)= inf sup E^[Z\ . (15) 

XeC~(s3,c2) yeC~(K3,c2) 

l!Z|li2(K3) = l. Z=(J) 

With the notation Z = (y) , we define the restricted min-max problem 

\i{v,B) inf sup E^[Z] . 

n-^x=o , o<||x||2 <i ^ ' \Y! 

^ i2(H3) „cy^o. I|V|I^2(h3)=1-|I^II^2(j.3) 

We show below that this restricted problem is actually a one-dimensional prob- 
lem. For this purpose, let us define the function a,^ : R x IR+ M given by 

Jo vs'^ + z'^ 



and implicitly define /i£ as the unique solution of 



l + fic[.f, i^, B]+va^[z) 



+ {l-ua^{z))\f{z)\'\dz 



\I{z)?dz 
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Theorem 3. For all B > and v £ (0, 1), 

Proof. The definition of Af (i^, S) is equivalent to 

Xf{iy,B)^ inf sup I ■(/;, 

t.T£fr;o' x.c-(M^c^) VV-*^3a.3 -1-^ 



or 



Af(z.,B) = inf sup 1 77^, 77 

Il"'/'lli2(.3,=l 



with the notation r = |a;| = ^ x\. For any given such that TT^c/) = 0, 
the supremum in x is achieved by the function 

Xc W\ = (Vi + (A£ J^, 5] + 1) Be) 'nd-i (T3 9.3) <^ , 
with Vc{x) -.^nc-^nc and 

1 - ^ - i 0-3 9^ 



A£[0, J/, B] sup 



77 V' 77-0 



Since (T3 = erf = I, this yields the expression 



*ec~(M3,c2) Jjj3 [Vi + A£[0, 1/, BJ + 1 V rJ 



dx 



nl4. = 0, ||0||_^2(e3, = 1 

Now, with the notations of Lemma |21 for all 7^ > 0, for all /i : M ^ M, we 
have /jj2 4>i h{r) (p*^, dxi dx2 = 0, with s = \/ x\ + Xj, r = ysM-aJf, and 

Af(z/,73) = inf inf ({Vc+Xc[cl>,i^,B] + l)-'d.,,cl>,d,,^)+(l-- 

ieN ii*ii_l2(k3)=i V V r 

i/>eCJJ° (M,span((/>^)) 

A simple calculation shows that for any g £ C^(M^, C^), 



£>0 

and also that 



r 



A simple integration by parts shows that for all i > 0, af < af_^ a.e. When 
minimizing, only the i = component has therefore to be taken into account. □ 
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Corollary 3. For all v £ (0, 1), the function [0, +oo) 3 B i-^ ^li^^; B) is nonin- 
creasing in B . 



Proof. A simple change of variables shows that 

aQ{z)— / — ds . 



— + z2 

g -t- ^ 

By Theorem |31 and according to the definition of /i£ [/, v, B] , this implies the 
monotonicity of Af (z^, •) in [0, oo). □ 

Proposition 10. For all B > 0, the function v i— *■ Af (i^, i?) is continuous in 
the interval (0,1) as long as it takes its values in (—1,1). Moreover, for any 
V E (0, 1), as long as \f{i',B) takes its values in (—1, 1), there exists a function 
Z G Range{n) with ||.^||l2(r3,c'') — 1 such that £[y[Z] = Af (z^, B). 

The proof is similar to the one of Proposition |2| □ 

The above proposition enables us to define 

BcH :=inf{B>0 : Xf {ly, B) = -1} . 

Recall that lim^^o+ '^+(^) = 1/2- We are now ready to state and prove the 
main result of this section. 



3.4- Asymptotic results. 

Theorem 4. Leti^e (0,t>). For any B e {l/d+{iy)^ ,mm {B{iy), Bdi^+i^^^'^)}) , 
we have 

Af fi^ + i^-'^/2 < Xi{ty,B) < Xf{iy-iy^/^,B) . 



Notice that the right hand side inequality holds for any B e (l/d+{i')'^ , B{iy)) . 
Proof. To prove the upper estimate, we use (I15f) and notice that 
Ai(t^,S)< inf sup S^Z), 

XeCg°(l3,c2) l'eC^(R3,c2) 

nlx=o,n^x^o ,x\ 

|Z|li2(K3)-l . ^-\y) 

since adding the condition W^X = increases the value of the infimum. Then, 
by Propositions IHl and 1^1 

Ai {u, B) < inf sup f ,_,3/2 [n Z] + [77= Z] 

XeCg°(M3,c2) Fgc„°°(R3,C2) 

£:,_,3/.[7TZ]-d_(;.)VS||77^r||2 



< inf^ ^ sup ^^.^^ ^iirTcvll2 



XeCo~(E3,c2) ^^^^-(R^C^) II^^IIl2(R3) + ll^£^llL2(R3) 

n'^x=o,n ^x^o u ^ / \ j \ j 

^ ■ r f,^_^3/2[7TZ] ^ , 

< mf sup II v"! , B) . 

n<j^x=a,n i^xi^a u / \ i 

Next, we establish the lower bound. By taking in H15() a smaller maximizing 
class of functions we decrease the maximum: 

Xx{v,B) > inf sup S^Z] . 

i7-y = 0. ||Z||^2(b3,=1 
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Therefore, by Propositfons|Hland|ni 

Ai [v, B) > inf sup ( £,+,3/2 [77 Z] + £,+^[n^ Z] ) , 

[n z] + d+{u) Vb \\n^cn%^,3^ 



nlY=a, ||z||^2(b3) = i 

Let us now notice that for every X e C^(M3, C^), X 7^ 0, 

is uniquely achieved at some Yc [X] because of the same concavity argument as 
in the proof of Theorem ^ after noticing that for any B £ {0,Bc{i' + i^^^^)), 
-1 < Xf{i^ + v^/'^,B) < \c[X, V + v^/^,B] and Vc > 0. Recall that 

Af(^/ + ^/'^/^B) = inf Xc[X,v + v^'^ , B] . 

x:ec^(K3,c2) 

Denoting Zc[X] = (y^x])^ given X, we find 



sup 



5,+,3/2 [n z] + d+{v) VB wn^cXWi.^^s, 



yeCo=^(R3,C2) ||77Z||^2(R3) + ||77^X||^2 



£,+,3/2(7/- Zc[X]) + d+iiy) VBWn^x 

- \\nzc[x]\\l2,^.,, + \\nix"^ 



2 



L2(K3) -r w-^J-c^^ llL2(R3) 

XcX, V + v^i\B\\\nzc\xw\2^^,^ + d^iv) Vb \\nix\\l^^^^ 
\\n Zc[x]\\\2(jf3) + ll-^£-'^lli2(K3) 



for B large enough so that d+{i>)\/B > Xf{h' + v^^"^ , B). As we shall see below, 
this is always possible. 

Note indeed that on (0, dj^{v) < ^2. Hence, d+{vy^ > 1/2. Now, by 
monotonicity (see Corollary OJ, Af (5, B) < Af (5, 1/2) for all S e (0, P + v^) = 
(0, 1 — v^/2) and for all B > 1/2. Moreover, one can prove very easily that 
Af((5,l/2) < 1 for all S e (0,1). Indeed, by Theorem|3 for aU B, X^{S,B) = 
inf^ Hc[f-, ^1 B]. A simple scaling argument shows that we can make J^\f'\'^dz 
as small as we wish while keeping J^\f\'^dz constant. Taking into account the 
definition of iic[f, 73], this shows that Af (^, 73) < 1 for all 73, for all 5 e (0, 1). 
Therefore, for all B > d+iv)-"^ , d+{v)^/~B > Af (i/ + v'^/'^,B) holds true. □ 

From Theorem 01 we deduce the following 
Corollary 4. Let v e {0,D). Then 

Bc{v + v^''^) < B{v) < Bc{v-v^''^). (16) 

From Theorem O and Corollary 0] better estimates of the critical magnetic 
strength B{v) than those of Theorem |21 can be established for v small. 
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Theorem 5. The critical strength B{v) satisfies: 

lim v\ogB{v) ~ TT . 

Proof. Because of Corollary 0] B(y) can be estimated using Bc{S) with S = 
V ± v^l"^ . This amounts to look for the smallest positive B for which 

A£(5,B):=1+ inf I {"^^^^wr^ ~ ^ a^{z)\f{z)\^ dz 

is such that 

Xc{S,B)^-l . 

Using the identity 




by the changes of variable and function 

y{z):= j\l{t)dt , f (^-fjj = B^/^ giy) , 
one tranforms the above minimization problem into 

WgU> l{]\9'iyr~S\giy)f) dy. 

Hence, 

\c{5,B) = 1 + Vb{\c{5,1)-1) . (17) 

For a given 5, let k — hi{5) := 5 [l — A£((5, 1)) and := l/aj(z(y)). The 

problem is reduced to look for the first eigenvalue Ei — Ei{5) of the operator 
—dy-\- k{S) fi{y), namely to find 6 such that 

The function aj satisfies 
aliz) < alio) 

There exists therefore a constant c > such that ii{y) < ce'^' for any y G R. 
To get an upper estimate of -Ei, we may now consider the function gi{y) := 
cos(7ry/2) on (—1, 1) B y and the rescaled functions ga{y) '■— <^^^^^gi{ ■ /(^)- 

EiiS)<^ + K £ |gi|V(ay)dy< ^ + Kc j\'^\y\\g,\^ dy 

< — + 2Kc(e'^-l) . 
Optimizing in a in the above expression, we choose <t = a{5) satisfying 

TT^ = 4 K C (T^ e"' , 

which implies 

cr((5) ^ — log K =: as ^ oo as (5^0. 



al{z) - — as |z| ^ oo 
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A Taylor expansion at next order shows that 

a{6) - as ^ -3 log {as) , 

which yields 

£^iW<^(i + o(i)). 

Next, in order to obtain a lower estimate of Ei{S), we consider the function 
fis which is equal to in the interval {—as, as) and equal to k ^{as) elsewhere. 
The function j^s is positive, but below the function KiJ,{y). Then observe that 
the first eigenvalue of the operator —dy + ^s{v), that we denote i?f , satisfies the 
equation 

^ias = arctan (\/^^^^^) , 
and as 5 goes to 0+ this implies 



E,{6)>Ei^^^{l + o{l)). 
Summarizing, what we have obtained is 

4(logK(,5)) 

So, imposing Ei{5) = (5^, we get k{5) = e~^(i+°(^'). Since by (fT7|l . 

-l = \^{6,Bc{5))=l + ^/B^)^, 

we get 

which, together with H16|l. concludes the proof. □ 



Appendix^: selfadjointness of Hb — | • | ^ 

With the notation P = —iW, consider a Dirac operator of the form 

Ho = a- P + mP + Voix) 

defined on (C(J°(R3))4. We assume for instance that {Yo)^,J{x) G Lf^^{^^)- If 
f €C^{R^), then the following identity holds on (C^(M3))4 . 

Hof - fHo-ia-V f = . (18) 

Still denote by Hq the a selfadjoint extension of a ■ P + m /3+Vo(x) with domain 
Vo and let / e C§°{R^) be such that < f{x) < 1, f{x) = 1 for |a;| < 1, f{x) = 
for \x\ > 2, and fa{x) := f{x/a). If 

V:^ {g = fa^P\a>l, ^ e Vo} , 

^ The following proof was explained to us by George Nenciu to whom we are grateful. 
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then in aU interesting cases, including the case of local Coulomb singularities 
with < 1 (see the characterisation of I?o in ll5 )) one has T) G Vq, and then 
by density, Identity (|18|l implies that for all G Vq, 

{Hofa- faHo-ia-V fa)^^0 . (19) 

A standard characterization of selfadjoint operators that we are going to use is 
the following: Let T be a closed symmetric operator. T is selfadjoint if and only 
if 

Ker(r* ± i) = {0} . 

The first remark is that Hq is essentially selfadjoint on V. Suppose indeed 
that 

(f , {Ho ± i) ifa^)) =0 Va>l, VVePo 
for some if e {L'^{R^))'^. Now from (O, 

{W, [Ho ± i) fai^) = (Ho ±i)^)- ((1 - /a)f , (Ho ±i)ij)+i {9, (a • VfaW , 
so that, taking a cx3 at fixed ip, one obtains 

(If, {Ho ±i)^)^0 V V G , 

and then tf' = 0. 

This result also applies to the case 

Vo = -^, 0<z/<l. 

See |13II17| for more details. Define ||| • ||| as the matrix norm. Suppose that Vi 
is locally L°° and more precisely satisfies 

M{R) := sup |||Vi(a;)||| < oo V i? G R+ , lim M{R) = oo . 

\x\<B. 

Consider on V the operator H = Ho + Vi . 

Lemma 3. Under the above assumptions, H is essentially selfadjoint on V. 

Proof. Let XR be the characteristic function of the set {x \ \x\ < 2R}. Since 
Xr^i is bounded, Hr = i?o + X_rVi is essentially selfadjoint on V. Suppose that 
there exists ^ G (L^(M^))'^ such that 

(<f , iHo± i) ifa^)) =0 Va>l, VtAgDo 

and assume that ||!f'||i2(j|3) — 1. Using (|18|l . one deduces that 

< fa^, {H± l) iP) = -I (If, a ■ Vfa^) . 

Observe now that < /a!f , {H± i) V^) = (/a!f , {Ha± i) V"), which amounts to 

(/af , {Ha± l) V-) - -I (f , a ■ Vfa^) • (20) 

Since Ha is essentially selfadjoint on V, Range(iJa± i) is dense in V and there 
exists ^/ij e 2? such that 

iHa±t)^Pi ^ fa^ + Sa (21) 

with ||(5a||L2(E3) < 1/a. Also notice that 

\\{Ha±i)^a\\h(R^) = \\Hai^a\\h{R^) + Hi\\h(R^)> II^?IIl2(R3) • 



Relativistic hydrogenic atoms in strong magnetic fields 



27 



From l)2()|l written for and (|21|l . we get 
SO that 

ll/a^lli2(R3) < ^ + l|V/a||L~(R3) 1 1 V'^ II (R3) 

< ^ + ||V/a||L-(R3) \\{Ha±i)lpi\\L-(^^) 

< ^ + l|V/a||L-(R3) + • 

For a CX3, ||/a'?'|lL2(R3) ^ 1 and faWh^^im?) 0, a contradiction: if' = 0. □ 
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